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Cavity opto-mechanics studies the coupling between a mechanical oscillator and a cavity field, 
with the aim to shed light on the border between classical and quantum physics. Here we report 
on a cavity opto-mechanical system in which a collective density excitation of a Bose-Einstein 
condensate is shown to serve as the mechanical oscillator coupled to the cavity field. We observe 
that a few photons inside the ultrahigh-finesse cavity trigger a strongly driven back-action dynamics, 
in quantitative agreement with a cavity opto-mechanical model. With this experiment we approach 
the strong coupling regime of cavity opto-mechanics, where a single excitation of the mechanical 
oscillator significantly influences the cavity field. The work opens up new directions to investigate 
mechanical oscillators in the quantum regime and quantum gases with non-local coupling. 
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Cavity opto-mechanics has played a vital role in the 
conceptual exploration of the boundaries between clas- 
sical and quantum- mechanical systems [1]. These fun- 
damental questions have recently found renewed interest 
through the experimental progress with micro-engineered 
mechanical oscillators. Indeed, the demonstration of 
laser cooling of the mechanical mode [H |31 IH ISl El IZ] 
has been a substantial step towards the quantum regime 

[SlISlIio]. 

In general, light affects the motional degrees of free- 
dom of a mechanical system through the radiation pres- 
sure force, which is caused by the exchange of momentum 
between light and matter. In cavity opto-mechanics the 
radiation pressure induced interaction between a single 
mode of an optical cavity and a mechanical oscillator is 
investigated. This interaction is mediated by the optical 
path length of the cavity which depends on the displace- 
ment of the mechanical oscillator. 

New possibilities for cavity opto-mechanics are now 
emerging in atomic physics by combining the tools of cav- 
ity quantum electrodynamics (QED) [llj il2j with those 
of ultracold gases. Placing an ensemble of atoms inside 
a high-finesse cavity dramatically enhances the atom- 
light interaction since the atoms collectively couple to 
the same light mode [El [HI [HI [16l [HI US] . In the dis- 
persive regime this promises an exceedingly large opto- 
mechanical coupling strength, tying the atomic motion 
to the evolution of the cavity field. Recently, a thermal 
gas prepared in a stack of nearly two-dimensional trap- 
ping potentials has been shown to couple to the cavity 
field by a collective center of mass mode leading to Kerr 
nonlinearity at low photon numbers |16j and back-action 
heating induced by quantum- force fiuctuations [19^ . 

A crucial goal for cavity opto-mechanical systems is 
the preparation of the mechanical oscillator in its ground 
state with no thermally activated excitations present, 
yet at the same time providing strong coupling to the 
light field. Here we use a Bose-Einstein condensate as 
the ground state of a mechanical oscillator and thereby 
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FIG. 1: (A) Cavity opto-mechanical model system. A me- 
chanical oscillator, here one of the cavity mirrors, is coupled 
via radiation pressure to the field of a cavity whose length 
depends on the oscillator displacement. (B) Coupling a Bose- 
Einstein condensate dispersively to the field of an optical 
high-finesse cavity constitutes an equivalent system. Here a 
collective density excitation of the condensate acts as the me- 
chanical oscillator which strongly couples to the cavity field. 
Feedback on the cavity field is accomplished by the depen- 
dence of the optical path length on the atomic density dis- 
tribution within the spatially periodic cavity mode structure. 
In contrast to opto-mechanical systems presented so far this 
mechanical oscillator is not based on the presence of an exter- 
nal harmonic potential (e.g. a spring). It is rather provided 
by kinetic evolution of the condensate density excitation. 



suppress thermal excitations of the oscillator to an un- 
precedented level. The cavity field couples to a collec- 
tive density excitation of the Bose-Einstein condensate 
which matches the cavity mode, resulting in an exceed- 
ingly large coupling strength. Despite the absence of an 
external restoring force for the mechanical oscillator, the 
framework of cavity opto-mechanics can be applied since 
only a single excitation mode of the Bose-Einstein con- 
densate is involved, see Fig. 1. 

In our experimental setup [HI [20] a Bose-Einstein con- 
densate of typically 1.2x10^ ^''Rb atoms in the |F, ttt.^) = 
|1, — 1) ground state is coupled to the field of an optical 



ultrahigh-finesse Fabry-Perot cavity. Our system is in 
the strong coupling regime of cavity QED, i.e. the maxi- 
mum coupling strength between a single atom and a sin- 
gle a~ polarized intracavity photon g^ — 2tix 10.9 MHz is 
larger than both the amplitude decay rate of the atomic 
excited state 7 = 27r x 3.0 MHz and that of the in- 
tracavity field K = 27r X 1.3 MHz. Trapping the con- 
densate within the cavity is accomplished by a crossed- 
beam dipole trap with trap frequencies {utxTi^yT^z) = 
2tt X (222,37, 210) Hz, where x denotes the cavity axis 
and z the vertical axis. The cavity has a length of 178 
pm and its TEMqo mode has a waist of 25 pm. The mode 
maximally overlaps with the condensate having Thomas- 
Fermi radii of {R^,Ry,R^) = (3.3,20.0, 3.5)iim. All ex- 
periments presented here were performed without active 
stabilization of the cavity length, which would give rise 
to an additional standing wave potential for the atoms 

[SlIIHlIig. 

The coupled dynamics of the Bose-Einstein conden- 
sate and the cavity field is driven by continuously ap- 
plying a weak pump laser field along the cavity axis 
(see Fig. 1). The light transmitted through the cavity 
is monitored using a single-photon counter and serves as 
a probe for the dynamics of the system. With a detuning 
of Aq = Ldp~Lda > 10'*7 between pump laser frequency ujp 
and atomic D2 line transition frequency uja spontaneous 
emission can be mostly neglected. 

Figure 2 shows the response of the system while scan- 
ning the pump frequency across the optical resonance. 
We observe a characteristic transmission signal (Fig. 2A) 
which exhibits a sharp rising edge and subsequently reg- 
ular and fully modulated oscillations (Fig. 2B) lasting for 
about 2.5 ms. These oscillations start at a frequency of 
about 37 kHz which slightly decreases over the train of 
oscillations and does not depend on the speed at which 
the pump frequency is varied. Similar responses of the 
system were measured for lower pump strengths at the 
same detuning (Fig. 2C) as well as for pump-atom de- 
tunings of up to Aq = 27r X 300 GHz, provided the pump 
rate was increased sufficiently. Moreover, when continu- 
ing the pump frequency scan we observe a second train 
of oscillations in the vicinity of the stronger coupling a~^ 
resonance [18]. This is in accordance with the observa- 
tion that the condensate remains intact during probing, 
which is directly inferred from absorption images taken 
subsequent to probing. The observed oscillatory behav- 
ior is obviously in strong contrast to a Lorentzian shaped 
resonance curve which would be expected for an atomic 
ensemble frozen inside the cavity, i.e. when the atomic 
external degree of freedom is neglected. 

To describe the driven BEC-cavity dynamics we con- 
sider a one-dimensional model in which the atomic mo- 
tion along the cavity axis is quantized. Justified by the 
large detuning between pump laser frequency and atomic 
resonance we adiabatically eliminate the internal state 
dynamics of the atoms. Denoting the creation opera- 
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FIG. 2: (A-C) Response of the continuously driven BEC- 
cavity system. Shown is a single trace of the cavity trans- 
mission (averaged over 2ps) while scanning the cavity-pump 
detuning at a rate of +2-k x 2.9MHz/ms across its cr^ reso- 
nance [18]. The pump rate corresponds to a mean intracav- 
ity photon number on resonance of 7.3 ± 1.8 (A and detail 
B) and 1.5 ± 0.4 (C). The photon count rate for one mean 
intracavity photon is 0.8 ± 0.2 MHz. The dead time of the 
single-photon counter is 50 ns which leads to a saturation of 
high photon count rates. The pump laser was blue detuned 
by Aa = 27r X 32 GHz with respect to the atomic resonance. 



tor for cavity photons by a' and the condensate wave 
function (normalized to the atom number N) by ■(/', the 
equations of motion for the coupled system read [HJ E?| 



ihil){x) 



^h^ <P 



(a}a)hUo cos^(fcx) 



2m dx-^ 



(1) 



ia = — I Ac — C/o(cos^(fca;)) + ZK ja + J77. (2) 

Here, T4xt denotes the weak external trapping potential 
for atoms with mass m, and giu the effective atom-atom 
interaction strength integrated along the transverse di- 
rections. 
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FIG. 3: Steady state and dynamical behavior of the BEC-cavity system in the two-mode model. (A) Mean intracavity photon 
number and corresponding oscillator displacement X versus the cavity-pump detuning Ac for the steady state solutions of Eqs. 3. 
The curves correspond to mean intracavity photon numbers on resonance of rf /n^ = 0.02, 0.07, 1 and 7.3, and a pump-atom 
detuning of An = 27r x 32 GHz. The inset highlights the bistable behavior for pump amplitudes larger than rj^r ~ 0.27k. (C-E) 
Evolution of the system in the mechanical phase space depicted for three subsequent situations (Ac = 27rx (200, 209.7, 215) MHz) 
corresponding to the markers in A and rf / k^ = 7.3. The stable and unstable steady state configurations are displayed as filled 
and open circles respectively. Dashed lines show representative evolutions for different starting conditions. Coloring indicates 
the modulus of the time evolution field {X,2P/h). The solid lines in D and E correspond to the experimental situation in 
Fig. 2A and show the evolution of the system while scanning Ac at a rate of 27r x 2.9 MHz/ms across the resonance with the 
system initially prepared in the lower stable solution. (B) Intracavity photon number |ap and corresponding transmission 
count rate r (including detection shot noise and averaging over 2 ps) for the system circling along the solid line in E. For 
integration of the equations of motion a coherent intracavity field a was assumed. 



Equation 1 describes the condensate dynamics in a 
dynamic lattice potential. Its depth is determined by 
the mean intracavity photon number {a''' a) which de- 
pends in a non-local and non-linear way on the con- 
densate wave function ip itself. For a single intracav- 
ity photon the potential depth is given by the light 
shift Uo = gg/Aa- The coupling between cavity field 
and atomic external degrees of freedom is mediated by 
the spatial overlap (cos^(fca;)) = / \tp{x)\'^ cos^{kx)dx be- 
tween atomic density and cavity mode structure, with 
wavelength A = 27r/fc = 780 nm. This mode overlap de- 
termines the effective refractive index of the condensate 
and with it the frequency shift of the empty cavity reso- 
nance in Eq. 2. The pump laser which coherently drives 
the cavity field at a rate 77 is detuned from the empty 
cavity frequency Wc by Ac = cjp — lUc- 

The observed BEC-cavity dynamics (Fig. 2) can be 
described in a homogeneous two-mode model where the 
macroscopically occupied zero-momentum state is cou- 
pled to the symmetric superposition of the ±2fifc mo- 
mentum states via absorption and stimulated emission 
of cavity photons. The corresponding wave function 



reads tp{x,t) — co{t) + C2{t)\/2cos{2kx) with probabil- 
ity amplitudes cq and C2 fulfilling |co(<)P + |c2(i)P = 
N. The mode overlap is then given by (cos^(fcx)) — 
(N + A/2Re(cQC2))/2. It oscillates under kinetic evo- 
lution of -0 at four times the recoil frequency Wrec — 
hk'^/{2m) — 27: X 3.8 kHz, with the atom-atom interac- 
tions being neglected at this stage. This leads to the 
natural definition of a harmonic oscillator with displace- 
ment X = 2-y/l/A^Re(coC2) in units of the oscillator 
length, and its conjugate variable P = h^J \ / Nlvii{d^C2) ■ 
The equations of motion (Eqs. 1 and 2) then read for 



|C2|V|C0P«1 



X + (4wi.c 



fx 

ia 



-{A + iK)a + ir] (3) 



and describe a mechanical oscillator coupled via the ra- 
diation pressure force to the field of a cavity whose reso- 
nance frequency shift A = A^ - UoN/2 - Uo/2y/N/2X 
depends linearly on the oscillator displacement X |23j . 
The coupling strength between optical and mechanical 
resonator can be varied via the atom-pump detuning Aq 



which allows us to experimentally enter the regime of 
strong coupling. 

From this equivalence to cavity opto-mechanics we can 
anticipate bistable behavior. Indeed, for pump rates 
larger than a critical value rjcr we find three steady state 
solutions for the oscillator displacement X, with two of 
them being stable (see Fig. 3A) [Ml Ell 1^- The sys- 
tem prepared below the resonance will follow the steady 
state branch until reaching the lower turning point, where 
a non-steady state dynamics is excited. This dynam- 
ics is governed by the time scale of the mechanical mo- 
tion since the cavity damping is two orders of magnitude 
faster. Thus we can assume that the cavity field follows 
the mechanical motion adiabatically and that retardation 
effects, underlying cooling and amplification, are negligi- 
ble [35]. Numerical integration of the coupled Eqs. 3 for 
our experimental parameters results in fully modulated 
oscillations of the cavity field and cavity output (Fig. SB), 
which is in very good agreement with the experimental 
observations (Fig. 2B). 

Further insight is gained by examining the dynamics in 
the phase space of the mechanical oscillator, spanned by 
X and P (Fig. 3C-E). Without cavity field the time evo- 
lution would simply correspond to a clockwise rotation 
at 4(j;roc- Yet, when photons enter the cavity the evo- 
lution is affected by light forces. This is the case along 
the vertical resonance line determined by the resonance 
condition A = 0, as shown in Fig. 3C-E (red line). 

Initially the condensed atoms are prepared at the sta- 
ble phase-space point {X, P) = 0, see Fig. 3C. Increasing 
the detuning A^ across the resonance renders the system 
instable and triggers parametrically excited oscillations, 
as indicated by the solid line in Fig. 3D. The evolution 
along this path is dominated by the free oscillator dynam- 
ics which gets periodically interrupted by the interaction 
with the cavity light field. Fig. 3D and E. This behav- 
ior is closely related to the matter-wave dynamics of a 
kicked rotor which is operated at an antiresonance where 
the accumulated phase factor between two kicks inhibits 
occupation of higher momentum modes |26j . 

The frequency of these oscillations decreases contin- 
uously over observation time (see Fig. 4). This is ex- 
pected when actively scanning the cavity-pump detuning 
Ac which shifts the resonance line in the phase space di- 
agram and leads to an adiabatic change of the system's 
circling path (compare Fig. 3D and E). 

A precise quantitative understanding of the observed 
frequency and its decrease is obtained when taking atom- 
atom interactions, the external trapping potential and 
atom losses into account. The atom-atom interactions 
result in a shift of the bare oscillation frequency 4a;rec = 
27r X 15.1kHz by the mean field energy, which in the 
Thomas-Fermi limit equals 4/7 times the chemical poten- 
tial /Lt = 27r X 2.4kHz [2H]. The trapping potential gives 
rise to a Fourier-limited broadening of the initial momen- 
tum distribution and accordingly introduces a damping 
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FIG. 4: Oscillation frequency while scanning over the reso- 
nance. The frequency within time bins of 50 ps was obtained 
from a peak-detection routine applied to the cavity transmis- 
sion data averaged over 10 ps. The data (filled circles) is an 
average over 23 traces referenced to the start of the oscilla- 
tions. The error bars indicate the standard deviation of the 
mean. Open circles show the result of a numerical integra- 
tion of the ID system taking atomic interactions and external 
trapping into account (Eqs. 1 and 2) [57]. The mean intra- 
cavity photon number on resonance was 3.6 ± 0.9. To fit the 
slope of the data the effect of a dynamically induced atom 
loss during the time of oscillations of 1.5 • 10^ /ms was added 
to the experimental frequency chirp of Ac — 2.9MHz/ms. 
The background rate of atom loss was measured to be 45/ms, 
and an atom number of (116 ± 18) x 10"^ was deduced from 
absorption images taken after the oscillations. 



of the free running oscillator dynamics. This suppresses 
a double peak structure in the transmitted light which 
would be expected at the onset of oscillations for the ho- 
mogeneous two- mode model (see Fig. 3D). An enhanced 
atom loss during the oscillations accelerates the observed 
frequency shift by a factor of 2. The numerical integra- 
tion of the full ID model (Eqs. 1 and 2) yields very good 
agreement with our data (Fig. 4). 

The quantitative agreement between experiment and 
semi-classical theory, together with the observation of 
very narrow peaks in the fully modulated cavity trans- 
mission, indicates that our system is well localized in the 
phase space of the mechanical oscillator. Using a sec- 
ond quantized picture where the Bose-Einstein conden- 
sate acts as the vacuum state of the mechanical oscillator 
mode, we have estimated the expectation value for ther- 
mal excitations in this mode. It is found to be below 0.01 
for a realistic condensate fraction of 90% [27 . This ex- 
tremely pure preparation of the ground state of a meso- 
scopic mechanical oscillator is possible since the cavity 
couples only to one specific excitation mode. Due to the 
high finesse of the cavity a single coherent mechanical ex- 
citation leads to a detectable shift of the cavity resonance 
by 0.7k. Entering this strongly-coupled quantum regime 
of cavity opto-mechanics promises to be ideal for testing 
fundamental questions of quantum mechanics [8] [9] [10] . 

From the perspective of quantum many-body physics 
we have investigated a Bose gas with weak local interac- 
tions subject to non-local interactions mediated by the 
cavity field. Experimentally it should also be possible 



to enter the strongly-correlated regime where local inter- 
actions dominate over the kinetic energy. In this case 
the non-local coupling is predicted to give rise to novel 
quantum phases [251 EHl 131] • 



Materials and Methods 

Quantum mechanical mapping to cavity 
opto- mechanics 

To obtain a fully quantized description of the coupled 
BEC-cavity system we start from its Hamiltonian in sec- 
ond quantized form. Its one dimensional version reads 
after elimination of the internal excited-state dynamics 
and in a frame rotating at the pump laser frequency LUp 

Eiiiia 



H = 



¥{x)(^- 



-h^ <f 



+ V^^t{x) 



2m dx^ 

+hUo coa'^{kx)d^aj i'{x) dx 

+Ha_a ~ ^.A^a'a — ihri{d — a* ) + H^. 

Here, VP^ denotes the creation operator of atoms with 
mass m, and a^ that of cavity photons with frequency 
LJc, wave vector k = 27r/A and mode function cos(fcx). 
The maximum light shift which an atom experiences in 
the cavity mode is given by Uq = g^/^a with the atom- 
photon coupling constant g^. The pump laser frequency 
is detuned from the empty cavity resonance frequency Uc 
and the atomic transition frequency Wo by Ac = cjp — Wc 
and Aq = bjp — LOa respectively. Decay of cavity photons 
at a rate k is accounted for by the term H^. 

In case of weak atom-atom interactions Ha-a and 
a shallow external trapping potential V^xt the BEC- 
cavity system can be mapped onto the generic Hamil- 
tonian of cavity opto-mechanics [32]. To this end we 
expand ^{x) into the two spatial modes (j)o{x) — 1 and 
4>2{x) — \/2cos(2fcx) which dominantly contribute to the 
BEC-cavity dynamics. The corresponding bosonic an- 
nihilation operators are denoted by Cq and C2- Apply- 
ing the Bogoliubov approximation cq = yN and taking 
(C2C2) ^ N into account we get 

H = 4/icji.ccC2C2 + /if — Ac + ,g(c2 + C2) Ja'a 

—itirj{d — a' ) + H,^. 

Due to the presence of the atoms we obtain a shifted 
cavity-pump detuning of Ac = Ac — ^UqN. The matter- 
wave mode 02 plays the role of a quantum-mechanical 
oscillator with its oscillation frequency 4ct;rcc, determined 
by the kinetic energy of this matter-wave mode. The os- 
cillator is coupled to the cavity field with a collectively 
enhanced coupling strength g = Uo/2^/N/2. To which 
extend quantum fluctuations play a role in the system 



is determined by the ratio g/K [35]. For our experimen- 
tal parameters we obtain g/n = 0.3 (0.6) for the (T~(cr+) 
transition. This coupling strength can be tuned via the 
pump-atom detuning A^ and the atom number N which 
allows us to experimentally enter the strong coupling 
regime of cavity opto-mechanics. 

From our coupling strength we can deduce the effec- 
tive oscillator mass rUeB in the cavity opto-mechanical 
model system [531 135]. There the radiation pressure cou- 
pling strength is given hy g = cUcO-i^o/L where L de- 
notes the length of the cavity and Oho = y/fi/i'^iTicS'^m) 
the harmonic oscillator length with oscillator frequency 
ujm = 4a;i.cc- With this we obtain an effective mass of 
nics = 0.01 ng. 



Ground state preparation 

We calculate the broadening of the overlap operator 
u = J dxcos^{kx)^^{x)^{x) caused by thermal depletion 
of the condensate. To estimate the number of thermal ex- 
citations in the mechanical oscillator mode we compare 
this broadening with the zero-point fluctuations of u for 
the mechanical oscillator being in its ground state. We 
restrict the discussion to the direction along the cavity 
axis since only density fluctuations along this axis con- 
tribute. The atomic annihilation operator ^(x) is split 
into a condensate part and a thermal part 

^{x) ^ ipo{x)^/N^ + y^ ipi{x)di. 

Here i^q denotes the condensate wave function obtained 
in the Thomas-Fermi approximation for Nq condensed 
atoms, and ipi the harmonic oscillator eigenfunctions in 
the external trapping potential. For simplicity we ne- 
glect the effect of atom-atom interactions on the ther- 
mal atoms. We evaluate the variance of u in a thermal 
state of temperature T where the uncondensed atoms 
Nt — N — Nq are distributed over the excited state levels 
according to a Bose distribution with chemical potential 
/i = 0. The total number of atoms N is kept fixed and 
the condensate fraction is given by Nq/N = 1 — {T/TcY 
with the critical temperature T^ calculated using the ex- 
ternal trapping frequencies. Beside autocorrelations in 
the density fluctuations we find for the variance Au^ 



Au^ = (u^) - {u9 



= 2N„Y.M^,{dld, 



i^O 



+ ^A4(aIa,)(a]a,) (4) 

where we introduced the matrix elements Mij = 
0.5{(pi\cos{2kx)\ipj). The first term on the right hand 
side of Eq. 4 originates from interference between ther- 
mal atoms and the condensate. For Nt < Nq this contri- 
bution by far exceeds the second term which corresponds 



to purely thermal density fluctuations at a wave vector 
2k. Comparing Au^ with the zero-point fluctuations of 
the overlap (0|(m - N/2)'^\0) = N/8 at T = gives an 
estimate for the number nx = ^^- of thermal excita- 
tions in the mechanical oscillator mode. Numerically we 
find a thermal occupation ut below 0.01 for a condensate 
fraction of 90%. 



Numerical integration of the semi-classical ID model 

The integration of the semi-classical Eqs. 1 and 2 in- 
cluding the external trapping potential and atom-atom 
interactions is done on a spatial grid of size 13A « 10 pm 
with 10 points per A/2. The cavity damping rate is two 
orders of magnitude larger than 4a;,oc which character- 
izes the time scale of atomic motion. Therefore the in- 
tracavity field dynamics is adiabatically eliminated and 
described by the coherent state variable a(t). Quantum 
fiuctuations of the cavity light field are not taken into 
account. The initial ground state solution for Ac be- 
low the optical resonance is obtained numerically using 
imaginary time propagation. Subsequently, this solution 
is propagated in real time with times steps of dt = 100 ns 
while scanning Ac across the resonance. 
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